Dynamic computed tomography manifestations of simulated wooden foreign bodies in blood-saline mixtures with variable concentrations and retention times

Diagnosing wooden foreign bodies (WFBs) using computed tomography (CT) is often missed, leading to adverse outcomes. This study aims to reduce misdiagnoses by exploring the density variation of blood-saline mixtures in ex vivo models. Twenty Cunninghamia lanceolata sticks, selected as WFB models, were randomly assigned to five groups: a control group (saline) and four experimental groups immersed in blood-saline mixtures with varying concentrations. The samples were then placed in a constant-temperature water bath at 36.8 °C. CT scans were performed in the lowest and highest density areas, and the volume of the low-density areas was measured at the post-processing workstation. Finally, the effects of time and concentration on imaging were analyzed, and fitting curves were generated. The blood-saline mixture concentration and time significantly affected the CT number in the three areas. WFB images changed dynamically over time, with two typical imaging signs: the bull's-eye sign on the short axis images and the tram line sign on the long axis images. Fitting curves of the CT number in the lowest density areas with different concentrations can quantify imaging changes. The CT number of the lowest density areas increased with time, following a logarithmic function type, while the CT number of the highest density areas exhibited a fast-rising platform type. The volume of the low-density areas decreased over time. The time of damage caused by WFBs and the influence of varying blood and tissue fluid contents at the damaged site should be considered in the diagnosis. Imaging changes from multiple CT scans at different times can aid in diagnosis.


Materials and methods
Study design. An ex vivo model was established, in which WFBs were immersed in BSMs of varying concentrations, simulating the internal environment. Using CT, we obtained the dynamic changes of the CT number for the lowest density area, the highest density area, and the volume of the low-density area of the WFBs at multiple predetermined time points.
A 12-year-old artificial Cunninghamid lanceolata fast-growing forest in the southwest of Hubei Province, with a diameter of approximately 15 cm, was selected and provided by the School of Forestry and Horticulture, Hubei University for Nationalities. Wood from the same annual ring area on the same layer of the xylem was cut into small wooden sticks measuring 40 mm × 2 mm × 2 mm. The wooden sticks were then dried at 50 °C to constant weight (absolute dry state). Density was calculated based on mass and volume, ρ≈0.38 × 10 3 kg/m3. Finally, all small sticks with statistically significant differences in length, width, height, and volume were excluded. The procedure's flowchart is shown in Fig. 1. www.nature.com/scientificreports/ Ex vivo blood saline mixtures model. A constant temperature water bath (Jintan Baita Xinbao Instrument Factory) was employed. An appropriate volume of water was placed in the bath, and the temperature was set to 36.8 °C. Plastic tubes with a 5 ml capacity and good sealing were filled with BSMs of different concentrations. The prepared wooden sticks were fixed using thin plastic sheets and suspended in the center of the tubes, respectively. Physiological saline served as the water component of the BSM. Five concentrations were established based on the proportion of whole blood contained in the BSMs, including the control group (saline group), trial (T) group T25 (containing 25% human whole blood), T50 group (50% human whole blood), T75 group (75% human whole blood), and T100 group (100% human whole blood), with four replicates for each group.
Scanning methods. Scanning equipment and parameters. A 64-slice spiral CT scanner was utilized (Philips Ingenuity 64 Slice CT Scanner, Philips, Netherlands). The tube voltage was set to 120 kV, tube current to 320 mA, rotation time to 2 s/cycle, display field of view to 25 mm × 25 mm, reconstruction layer thickness to 1.0 mm, pitch to 0.984, and a standard algorithm (Filtered Back Projection) was applied for reconstruction.
Scanning. After establishing the test tube models, they were placed on the CT scanning bed for the initial CT scan to acquire baseline image information. Then, the models were submerged in test tubes containing BSMs of varying concentrations and placed in a constant temperature water bath. At 6-h intervals, each group of test tubes was removed from the water bath, scanned horizontally and vertically to collect data at each time point, and scanned up to 612 h (25.5 days) after model establishment 17 .
Image data acquisition and analysis. The raw data were uploaded to the GE ADW4.7 workstation for postprocessing analysis. Regions of interest in the highest density areas and lowest density central areas were outlined separately to obtain CT numbers. The workstation's volume measurement software was used to measure the volume of the WFBs' low-density areas, defined as an area with a CT number < 0 HU. After completing data collection for the four group replicates, the mean was calculated and entered into the database for analysis.
Statistical analysis. Repeated-measures analysis of variance and least significant difference (LSD) were employed for statistical analysis to explore the differences between the CT numbers of the lowest density areas, the highest density areas, and the volume of the low-density areas of WFBs at different times and BSMs with varying concentrations (P = 0.05). An aspherical test was performed on the variables. The differences between BSMs with different concentrations at the same time were tested for significance using a one-way analysis of variance and LSD (P = 0.05). All data were subjected to the Shapiro-Wilk normality test before analysis, and data that did not conform to normal distribution were transformed based on log (x + k) or 1/x. SPSS 25.0 software was used for data analysis and fitting, while Origin 2021 software was utilized for graphing.

Results
Dynamic image signs of WFBs changing over time. In this study, we observed two typical imaging signs by analyzing the CT scan images of WFBs. The first sign was the bull's-eye sign. In the short-axis images, the edges of the WFBs displayed ring-shaped high-density shadows, while the central areas showed roundshaped low-density shadows. Against the background of the high-density areas at the edges, they resembled bull's-eyes (Fig. 2). The second sign was the tram line sign. In the long-axis images of the multi-planar reconstruction, the WFBs exhibited strip-shaped high-density shadows on both edges and low-density shadows in the center, resembling a track (Fig. 2).

Effects of BSMs with different concentrations on CT number of the lowest density areas of
WFBs. Both BSM concentration and time had a highly significant effect on the CT number of the lowest density areas of wooden foreign bodies (P < 0.001, Table 1), with a highly significant interaction between them ( Table 1, P < 0.001). Overall, the increase in CT number in the lowest density areas decreased with increasing BSM concentrations in unit time. The CT number curves of each group of time-lowest density exhibited a rapid rise initially, followed by a slow rise (Fig. 3). Before the experiment began, there was no difference in the CT number of the lowest density areas of the wooden sticks between the groups. After the sticks were immersed in BSMs with different concentrations, the CT number of the lowest density areas of each group increased significantly within the first 6 h, and the CT number of the lowest density areas between the groups showed differences. The CT number of the low-density areas in the Saline group increased the fastest, 589 HU higher than the initial value (− 861.5 HU). The CT number of the lowest density areas in the T100 group increased the slowest, only increasing by 407.5 HU compared with the initial value (− 837.5 HU). Afterward, the changes in CT number in the lowest density areas of each group displayed a slowly increasing trend. Among them, the lowest density areas of the NC group and T25 group disappeared at the 90th and 234th hours, respectively (CT number of WFBs > 0 HU). The fitting curves of different CT numbers of each group over time are shown in Table 2.

Effects of BSMs of different concentrations on CT number of the highest density areas of
WFBs. Both BSM concentration and time significantly affected the CT number of the highest density areas (Table 1, P < 0.001). During the experiment, the CT number curves of the time-highest density areas exhibited a fast-rising platform type (Fig. 4). Before the experiment started, the WFBs in each group were all low-density without high-density areas; after being immersed in BSMs, high-density areas appeared at the 6th hour, and the Table 1. Time, concentration, and the interaction of the two on the CT number of the lowest density area, the highest density area, and the volume of the low density area repeated measurement analysis of variance results.

Effects of BSMs of different concentrations on the volume of the low-density areas of
WFBs. Both BSM concentration and time had highly significant effects on the volume of the low-density areas of WFBs (Table 1, P < 0.001), with highly significant interactions between them (Table 1, P < 0.001). The degree of the decrease in the volume of the low-density areas decreased as the concentrations of the BSM increased. The time and low-density areas volume curves of the Saline group and T25 group exhibited a rapid decline, while the time and low-density areas volume curves of the T50 group, T75 group, and T100 group decreased slowly (Fig. 5). Initially, there was no difference in the volume of the low-density areas among the five groups. After being immersed in the BSMs, the volume of the low-density areas of each group decreased significantly within the first 6 h. The degree of volume reduction in the low-density areas gradually decreased from the Saline group, T25 group, T50 group, T75 group, to the T100 group. Compared with the initial value (167.85 mm3), the volume of the low-density areas in the Saline group was reduced by 42.85 mm3; the volume of the low-density areas in the T100 group was only reduced by 25.09 mm3 from the initial value (167.59 mm3). Among them, the volume of the low-density areas disappeared in the Saline group at the 96th hour and the T25 group at the 240th hour (CT number of WFBs > 0 HU). Ultimately, the remaining groups still had low-density areas of varying sizes.

Discussion
Principal findings. This study demonstrates that WFB images change dynamically over time, exhibiting two typical imaging signs: the bull's-eye sign on short-axis images and the tram line sign on long-axis images. We developed fitting curves of CT numbers in the lowest density areas across different concentrations to quantify the imaging changes of WFBs. Both the BSM concentration and time have a significant effect on the CT number   Potential mechanisms. The effect of the wood's water absorption capacity on density changes. WFB density, which determines the CT number, is closely related to the water absorption capacity of wood materials 18 . This capacity primarily depends on the cellulose and pore structure of wood, as confirmed by Gao 19 . Within the first 6 h of this study, the CT number of the lowest and highest density areas of each group increased rapidly. This increase can be attributed to the large internal cell cavity pores of the wood during the initial stage of water absorption, which rely on the capillary system composed of axial tracheids to absorb water with the help of capillary force, while lignocellulose quickly absorbs water 20,21 . Moreover, it has been reported that the isotherm adsorption process of wood mainly involves rapid monolayer adsorption on the adsorption sites within the wood. This process may represent another significant reason for the swift density changes of WFBs within the first 6 h 21 . As relative humidity increases, second-layer and multi-layer adsorption gradually appear, resulting in more water molecules being adsorbed by WFBs and a corresponding increase in density 22,23 . As the water content of lignocellulose rises, its hydrophobicity also increases, causing the water absorption capacity to slow down 20 . The water absorption process of WFBs progresses gradually from the surface to the center; the further from the wood surface, the slower the water absorption rate, and the relatively slower the density change. Consequently, the density of the area near the wood surface of the WFB changes most rapidly, whereas the density of the central area changes more slowly (Figs. 1, 2).

Effects of BSMs of different concentrations on the density changes of WFBs.
The change in WFB density essentially represents a process in which the air within the wood gap, under the absolute dry state of the WFB, is continuously replaced by free water and other substances from the surrounding environment. Due to the different concentrations of BSMs, free water enters WFBs at varying speeds, which delays the rate of density changes 24,25 . Simultaneously, as BSM concentration increases, so does the viscosity of the medium, hindering the diffusion of water molecules and reducing the speed at which various molecules from the WFBs' surrounding environment enter the interior, thus delaying the density change rate of the WFBs 26,27 . Additionally, the complex composition of BSMs means that organic macromolecules and various cells in the blood may adhere to the surface, forming a barrier film that slows down the entry of water molecules into the WFBs 28 . Another study has indicated that the gas within the wood foreign matter and the water molecules entering it form a gas-liquid plane, which is not conducive to the direct exchange of water molecules and gas 29 .
Image characteristics of changes in the density of WFBs and for best practice. As the edge of the WFB is in direct contact with the surrounding environment, the air inside it is replaced at a significantly higher speed than the interior. Consequently, high-density areas appear at the edges while the interior remains low-density. The tram line sign is generally observed on the long-axis of the planar reconstructed image, as seen in the WFBs cases reported by Mertel et al. 30 . The bull's-eye sign is mostly visible on short-axis images, as demonstrated in the cases reported by Peterson et al. 1 . When imaging a WFB, it is essential to consider the patient's description and the position within the image to make an accurate diagnosis.
In this study, the CT number of WFBs in the body varies considerably over time. The initial CT number of WFBs upon entering the body ranges from − 800 HU to − 900 HU, with a density close to that of gas, necessitating differentiation from gas accumulation within the body. As WFBs remain in the body for longer periods, their edge www.nature.com/scientificreports/ density increases while the central area remains low-density, and the CT number of the low central area gradually increases until the WFB becomes high-density overall 30 . The CT number of the WFB progressively rises to − 100 HU to − 200 HU, at which point the density is close to that of fat, requiring differentiation from fatty tissue. After some time, the density of WFBs will be similar to that of water, soft tissue, blood, or calcification, necessitating differentiation from these counterparts [31][32][33][34][35] . WFBs typically exhibit rapid changes at the edge and slower changes in the central area, resulting in a ring-shaped high-density image. This feature can be used as a characteristic manifestation. Radiologists can calculate the CT number of wooden foreign bodies in different tissues based on the fitting curve equations of time and various concentrations to reduce the likelihood of missed diagnoses.

Strengths and Limitations.
This study presents a novel approach for investigating WFBs by considering the dry state as the initial state and dynamically observing the imaging manifestations and change patterns of WFBs. Additionally, we discussed potential factors influencing the imaging changes of WFBs and elucidated the inherent nature of these changes. This work contributes to a better understanding of WFB imaging manifestations and alterations. Unlike previous studies on WFBs, which were conducted either in vitro or in animals and examined at a single time point 6,7,35 , our research offers a more comprehensive perspective on WFB imaging. This is particularly important given the continuous changes in density and imaging appearance of WFBs in the body. However, this study has several limitations. First, we only investigated Chinese fir materials, and the CT number changes might differ for various types of wood foreign bodies. Second, we focused on a low-density wood, and it remains to be determined whether high-density wood foreign bodies exhibit similar patterns. Finally, the composition of the organizational microenvironment in the human body is more complex than BSM environment. Future research should aim to establish models to better replicate the in vivo environment. Despite these limitations, our study represents a significant advancement in the field and provides valuable insights for future investigations.

Conclusion
In summary, this study elucidates the influence of blood concentration and retention time on the density of WFBs within the BSMs. The findings provide valuable insights for the identification and differentiation of WFBs in clinical scenarios. Accurate diagnosis of WFBs necessitates a comprehensive evaluation of the patient's medical history, the surrounding BSM environment, image positioning, and temporal changes.
For cases where WFB retention is suspected, a follow-up scan after a certain period may be recommended. Characteristic alterations in the images from the two scans can facilitate more precise clinical diagnoses. Future investigations should focus on high-quality, well-designed studies and trials that compare multiple diagnostic tools to further enhance the understanding and management of WFB-related cases in clinical practice.

Data availability
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.